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We analyze the B — > cj>K* polarization puzzle in the Minimal Supersymmetric Standard Model 
(MSSM) including the neutral Higgs boson (NHB) contributions. To calculate the non-factorizable 
contributions to hadronic matrix elements of operators, we have used the QCD factorization frame- 
work to the Qs order. It is shown that the recent experimental results of the polarization fractions 
in _B ^ (f)K* decays, which are difficult to be explained in SM, could be explained in MSSM if there 
are flavor non-diagonal squark mass matrix elements of 2nd and 3rd generations, which also sat- 
isfy all relevant constraints from known experiments [B Xs"f,Bs n'^fi^yB Xsfi^ , B — ^ 
XsQ, AMs, etc.). We have shown in details that the experimental results can be accommodated with 
the flavor non-diagonal mass insertion of chirality RL, RL+LR, RR, or LL-I- RR when the NHB 
contributions as well as 0{as) corrections of hadronic matrix elements of operators are included. 
However the branching ratios for the decay are smaller than the experimental measurements. 



I. INTRODUCTION 



The recent experimental results for polarization fractions in 
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for the mode (j)K*'^. The amphtudes |^o| and | are longitudinal and transverse polarizations of decay amplitudes 

in the transversity basis which satisfy 

E i^^i'^i- (3) 

j=0,_L,|| 

The results deviate significantly from the SM prediction 

|AoP^1-0(1/to,), (4) 

n 

based on the naive counting rules which follow from a helicity argument ^ . This significant deviation is referred as a 
puzzle or anomaly in the literature. It has attracted many interests in searching for possible theoretical explanation 



in SM and new models beyond SM 



The naive counting rules are obtained with the naive factorization in calculating hadronic matrix elements. It may 
be possible to explain the data if including the as corrections to hadronic matrix elements in SM. It is shown that 
one can obtain \ Aq\'^ ~ 0.5 due to the annihilation enhancement from the integral containing end-point singularity 
in QCDF approach 7]. However, it is at the issue that makes the approach less-predictable. Li and Mishima point 
out that annihilation contributions are not enough to make \Aq\'^ ^ 0.5 H in PQCD factorization approach 0|. 
The effects from the final state interaction (FSI) have been studied in refs. |lo[[ll|. One can get \Ao\^ ^ 0.5, but 
|^oP(^ ~^ P^* \ < l^oP(5 ^ (f'K*) 10] which does not agree with the measurements . Moreover, it has been 

shown in ref. [ll| that such FSI effects would lead to \Aq\'^ : \A\\ p : \Aj_\'^ = 0.43 : 0.54 : 0.03 which clearly contradicts 
the data. Therefore, one may draw the conclusion that it is difficult to explain the data within the SM. 

A lot of works have been done to investigate polarizations of B ^ (pK* in models beyond SM. A model with 
right currents can give |Aop ~ 0.5 but simultaneously leads to ^ which is not in agreement with the 

data ya It is also shown that the RL or LR-I-RL insertion in MSSM can lead to \Aq\'^ ~ 0.5 due to the Csg 
enhancement, compared with that in SM [ij . However, wrong formulas for the as order hadronic matrix elements 
of the chromomagnetic dipole operator Qsg in the case of transverse polarization are used in Ref. jisl ] . As shown in 
refs. p. the as order hadronic matrix elements of Qgg for transverse polarizations are very small. Moreover, the 
neutral Higgs boson (NHB) contributions are not considered in the workjisj. Yang et al. show that the R-parity 
violating SUSY might explain the puzzle 20]. A model-independent analysis for contributions of new operators, i.e., 
the operators beyond the operator basis in SM, has been carried out in ref. Recently, an analysis of polarizations 
in the model with scalar interaction of tree-level flavor changing neutral current (e.g., the model HI two Higgs doublet 
model) has also been performed 21]. In this paper, we shall perform a detailed analysis of polarizations in _B — > (pK* 
as well as the decay rates in MSSM including neutral Higgs boson contributions and the as corrections of hadronic 
matrix elements. 

For the — > s transition, besides the SM contribution, there are mainly two new contributions arising from the 
strong penguins and neutral Higgs boson (NHB) penguins with the gluino and squark propagating in the loop in 
MSSM. The former is not important because the Wilson coefficients of QCD penguin operators in MSSM are not 
changed significantly, compared with those in SM. Although Csg can get a significant enhancement 

00, the 

hadronic matrix elements of Q^g in the case of transverse polarization are very small. The latter induces scalar 
operators as well as tensor operators due to renormalization. As well known, the effects of these new o per ators 
to leptonic Bs decays are significant [2^ , and their effects to some hadronic B decays are also important [2^ . 
For B — !■ VV decays, it is expected that the hadronic matrix elements of scalar and tensor operators can enhance 
transverse polarization fractions. Moreover, although the effects of the primed counterparts of the usual operators are 
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suppressed by nis/mi, and consequently negligible in SM, their effects in MSSM can be significant because they have 
the opposite chirality and the flavor non-diagonal squark mass matrix elements are free parameters which are only 
subjective to constraints from experiments. In particular, as discussed in ref. [isl . the primed counterparts of the usual 
operators have contributions to longitudinal and transverse polarizations different from those of usual operators and 
consequently could enhance the transverse polarization fractions. The relevant Wilson coefhcients at the mw scale 
have been calculated by using the vertex mixing method in Ref. [2^ and the mass insertion approximation (MIA) 
method in ref. ^fil]. In this paper we shall use the results given in ref. [2^ . For the hadronic matrix elements of 
operators relevant to the decays B VV, we shall use the BBNS's approach (QCDF) to calculate the ag order 
corrections to the naive factorization results. 

We show that polarization fractions of the decays can agree with experimental data within la deviation in MSSM 
with the parameter space satisfying all the constraints from Bs — Bg mixing , B Xg'j, B Xsg, B Xs/i+/i^ and 
Bg fi^ . In particular, the puzzle for polarization in _B ^ (j)K* can be explained, while not in contradiction to 

the measurements of other two vector final states, in quite a large region of parameter space because we have included 
the contributions of the primed counterparts of usual operators and NHB induced operators in MSSM with the as 
corrections of their hadronic matrix elements included. 

The paper is organized as follows. In Sec. II, we give the effective Hamiltonian responsible for the b ^ s transition 
in MSSM. In Sec. Ill, we present the decay amplitudes. In particular, the hadronic matrix elements of NHB induced 
operators to the order are calculated. The Sec. IV is devoted to numerical results. We draw conclusions and 
discussions in Sec. V. 

II. EFFECTIVE HAMILTONIAN 



The effective Hamiltonian for 6 — > s transition can be expressed as[2^ |2; 

^ p—u,c ^ 2— 3,..., 16 

+C77 Q77 + C'sg Qsg + C77 Q'tj + C'sff Q'sg ^ + li-C- (5) 
Here Qi are quark and gluon operators and are given by * 

Ql = {SaPl3)v-A{Pl3ba)v-A, Q2 = (SqPq) V-A(P/3&/3)y-A, 

Q3{5) = {Saba)v-A^^{qf3qf3)v-(+)A, Q4{6) = iSabf3)v-A^^{qpqa)v-{ + )A, 

1 q 

Q7{9) = l^{Saba)v -A'^eq{qpqi3)v + {-)A, Qs.[W) = 2^Sahi3)v-A'^eq{qj}qa)v+(-)A^ 

q 9 



Qii(i3) = (s^)s+py] —{qq)s~(+)p 1 

9 

<5i2(i4) = (si&j)s+P^ —{q^qi)s-{+)p-, 
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* For the operators in SM we use the conventions in Ref. l28l where Q\ and Q2 are exchanged each other with respect to the convention 
in most of papers. 



Qi6 = Si cr'""{l + 75) bj a^,y{l + 75) Qi , 

Q77 = -z-^rnbSaCT^" F^^{1 +75)&/3, 

= {^^bScCr^^^'G^^il + 75)&;3, (6) 

where {qiq2)v±A = 9i7''(l ± 75)^2, {qiq2)s±p = ± 75)92 \ p ^ u,c, q ^ u,d,s,c,b, Cq is the electric charge 
number of q quark, is the color SU(3) Gell-Mann matrix, a and (3 are color indices, and F^^, {G^^) are the photon 
(gluon) fields strength. The primed operators, the counterpart of the unprimed operators, are obtained by replacing 
the chirality in the corresponding unprimed operators with opposite ones. 

For the processes we are interested in this paper, the Wilson coefficients should be run down to the scale of 0{mb). 
Ci — Cio are expanded to 0{as) and NLO renormalization group equations (RGBs) should be used. However for the 
Cro and C7-Y, LO results should be sufficient. The details of the running of these Wilson coefficients can be found in 

^ nn 

Ref. 29]. The one loop anomalous dimension matrices of the NHB induced operators can be found in refs. \2(i l3(l| . 
There is the mixing of the new operators induced by NHBs with the operators in SM. The leading order anomalous 
dimensions have been given in 

Refs. 00 

The mixing of NHB induced operators with the chromo-magnetic operator 
can enhance the Wilson coefficient Cgg significantly Is^ l . Because at present no NLO Wilson coefficients cl'"* , 
i=ll,...,16, are available, we use the LO running of them in this paper. 

III. THE DECAY AMPLITUDE AND POLARIZATION 



We use the BBNS approach |3,|23| to calculate the hadronic matrix elements of operators. In the BBNS approach, 
the hadronic matrix element of an operator in the heavy quark limit can be written as 

{ViV2\Q\B) = {ViV2\Q\B)f [l + ^r„<] , (7) 

where (V1V2IQI-B)/ indicates the naive factorization result. The second term in the square bracket indicates higher 
order as corrections to the matrix elements We calculate the hadronic matrix elements to the as order in this 
paper. In order to see explicitly the effects of new operators in the MSSM, we divide the decay amplitude into three 
parts. The first one. Ho, has the same form as that in SM, the second. Ho', is for primed counterparts of the SM 
operators, and the third, Hn, is new which comes from the contributions of Higgs penguin induced operators. That 
is, we can write the decay amplitude for B — > VV as 

A{B ~^ViV2) = ^H^ 
v2 

= H^ + H^, + H^ . (8) 
The helicity amplitudes can be obtained by set A = 0, +1, —1 in above expressions, respectively. 



t Strictly speaking, the sum over q in expressions of Qi {i=ll,...,16) should be separated into two parts: one is for q=u, c, i.e., upper type 
quarks, the other for q=d, s, b, i.e., down type quarks, because the couplings of upper type quarks to NHBs are different from those 
of down type quarks. In the case of large tan/3 the former is suppressed by tan~^ 13 with respect to the latter and consequently can be 
neglected. Hereafter we use, e.g., C^^ to denote the Wilson coefficient of the operator Qii = {Sfe)s+p j^{cc)g_p. 
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A. Helicity amplitude 

Let e'^''' (i=l, 2) be the polarization vector of vector meson Vi, Ai = A2 = A in _B — > V1V2 due to the angular 
momentum conservation. The helicity amplitudes (A = 0, +1, —1) of i? ^ (j)K*'^ in MSSM are given by 



A 



H 



1 



2 



'A, 



(9) 
(10) 

(11) 



^T(i+75)(A) - -(at2 +r2a'i^2)Av'+A)(A) [>(12) 



where 



A 



T(l-75) 



r2 



(V+A) 



Av±A)(A) = (</>(e(^\g)|57^(l±75)s|0>(i^*(e('\pKO|S7^(l±75)fe|S(Pi3)>. (13) 

In ea. lfTnjl the coefficients af, i = 3, 4, 5, 7, 9, 10 have been given in ref. 0,0,0,0- Because there are contradicting 
results on penguin insertion contributions, especially the Csg effect to transversely polarized amplitudes, we revisit 
this part and confirm the results in ref. 0, Q| . We calculate coefficients in eq. H12(l and results are 



n2 



ai4 



'■15 



— {C12 
rrib 

rrib 



nrib 

TO, 
TOf, 



-{Ci 
{C 
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Cii , tts Cp Cii 



^[1 



47r 



■(^3+^^13)]} + 



47r 



C16 Ci4 OsCf , I M I Q^'' ^16 pA 



15 



[1- 



47r 



(^15 + ffl^s)]} + 



«s Ci5 

47r 



(14) 



where 



pA 
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pA 
J 1 ■ 



^(^ln^-G^(0)) + lln!^-G^(l) 



^2, 



-8 



-21n."''^ ^^A, 



G^A — GFm (1) 



'13' 



A' = A^3 + 4G, 



pA _ 
1 fi — 



-13 



pA 



for A = 0, 



li^m^ for A = ±l. 



(15) 
(16) 



In eq. p 5116(1 . we have defined 
Av-a{\ = 0) 



^T(l+75)(A = 0) 



G^(,s)= / da;G(s-ze,l-x)$^i(x) 
"'0 

Jo ^ Jo ^ 



G{s,x) =-4: dtt{l-t)h-i[s-t{l~ t)x] 



/ dx GF{s-ie,x), GF{s,x) = 

Jo ^ Jo 



dt\n[s-xtt\ (17) 
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with X = 1 — X. Here the distribution amphtudes of 4> meson are given by 



(18) 



where (f>\\ , g^'' , and 



$= / dyi<t>\iiy)~g^f\y)) (19) 

are defined in ref . ,3J| . Using the Wandzura-Wilczek approximation [s^ , one has 

G^±=0. (20) 

And numericaUy Gi^^(l) is smaller than GF^{1) by about a factor of two and GF^{1) = 0. Thus, the penguin 
contract contributions of (5i=i3, 15,16 to transverse amplitudes are smaller than those to the longitudinal amplitude. 
However, the penguin contract contribution of Qn to transverse amplitudes can be larger than that to the longitudinal 
amplitude, as Gt(l) is larger than G^{1). 



In eq. (|14(l . V"^ and Hj^, 1 are vertex and hard-spectator scattering contributions respectively and numerically not 



important. 

The amplitudes in transversal basis [35j for B VV are related to the helicity amplitudes by 

H- +H+ _ H- -H 

And the longitudinal polarization is defined by 



Ao=HO, A,|=i^^, A, = f^^. (21) 



h - . , ,7;. , (22) 



Using the identity 



we have 



B. Form factors for B <j>K* 



o^'^-i^ = -'-e'^'P'^ap,, (23) 



(0|s-a'-(l±75)s|O) = (5^V"T^e^''''")(0|5ap.s|O), 
{K*{e^^\pK')\sa^^''{l±l,)h\B{pB)) = (s^V" T ^e^'''"^)(if*(e(^',PK.)k-a,.6|i3(ps)). (24) 

Defining 

{K*{e(^\pK')\-sa^,h\B{pB)) - -ie^,p,e*''(p^Gi(s) + p^.G2(s)), (25) 
where s — and q= ps — Pk* , one has 

(X*(e(^),PK.)|5(T''"'Z.^|S(ps)) - ie^''""'e,,pBpPK'ACi + C2). (26) 
Comparing with the usual definition, one has (Gi + G2) = 2T'i. From cas. H24|l and p5|l . it is easy to obtain 
{K*{e^^\pK>)\-sa>^^q,{l ± ^^)b\B{pB)) - ze'''''"^e:pBpP/f*.G+ ± \e*^ [(m| - m2,.)G+ + sG_] 

T^e*-Ps(p^G++<z^G_), (27) 



where C± = Ci ± C2 , s = , and p = ps + Pk' ■ From eq. (|27|l , it follows that there are only two independent form 
factors in the matrix element of the tensor operator between pseudo-scalar and vector meson states. That is, we need 
not introduce three form factors in the matrix element, as done in the usual definition in ref. [4^ . Comparing with 
the usual definition of the same matrix element, one has 



1 



C+ ^ C_ = 273, T, = -{C++ . C- ) . (28) 



Define 

{0\sa^''s\<l>{e(^\q)) = f^ie^^^%- ~ e^^^" q^^), (29) 

we have the naive factorization amplitude of tensor operators as follows. 

^r(i±^,)(A) ^ q)\sa,A'^±^5)m{K*(.e^^\pK')\sa^-(.l±l5)b\BipB)) 
= 4 /^^ e^,p„e*/erpPsp],, 2T, ± (ej • e^) (m| - mj,, ) 



T{el ■ Pb) [el ■ pb) Ut^ + -j^^^T^) ] , (30) 
where the superscript (A) has been suppressed in the right hand of eq. (|30|l . Therefore, it follows that 

-4t(i±^,)(A = 0) = ^2^^\{ml-ml-ml,)T2- Am%pl ( + , ""^ , Tg 
-4t(i±^,)(A = +1) = T^f^ {{ml-ml.)T2T2mt,PcT,}, 

-4t(i±75)(^ = -1) = T^fi {{m%-ml.)T2±2mbPcTi} , (31) 

with pc is the center mass momentum in the B rest frame. 

The decay constants and the form factors of vector and pseudoscalar mesons are defined as usual 

(O|s7M(l±75)s|0(g, e(^))) = tUm^e^^\ (32) 

vriB + mfc* 
TUB + mK* 

±iq^.e(^> ■Pb'^^{A^{s)-A^{s)), (33) 
s 

where p = + pk* and q = pB — Pk* ■ The above equations lead to 

Av±a{X) = ef^)\-s^^{l±^^)s\0){K*{pK',e^^'^)\-s-i,{lTimB{pB)) 



(A)*M (A)*i/ p c 2y(m^) |.^), \ A ( 2\ 

/.i/ptrej ' PbPk' ±4 ■ e\ ' {niB + mK')Ai{m^) 



rriB + rriK' 



(A)* (A)* M^l) 

■PBe-2 -Pb \ 

mB + rriK' 



^From ea. (|34|l . we have 

-4y±A(A = 0) = Ti^^^ - vi]^. - m^^){mB +mK') Ai- 'im%pl 

-4v±a(A = +1) = ^f^m^ <{mB + niK') Ai^2'mBP. 



mB + rriK' 

V 



(34) 



niB + ran* 

^y±^(A = -1) = ^f^m^ {{tub + ruK') Ai±2mBPc 7- \ (35) 
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Comparing eq. (|31|l and eq. H35|l . one has 

Av^A ^ -^T(l±75) ■ (36) 

That is, the contributions of tensor operator are enhanced by a factor of mB/m^^ compared with those of vector-axial 
vector operators. Therefore, the contributions of NHB are sizable although there is a suppression factor mg/m\, in 
eq. (jH) . 



IV. NUMERICAL RESULTS 



A. Constraints from experiments 



We impose two important constraints from B — > Xg^ and Bg — > fi'^fi^ . Considering the theoretical uncertainties, 
we take 2.0 x 10^'' < Br(_B — > Xg'^) < 4.5 x 10^'', as generally adopted in the literature. Phenomenologically, 
Br(i3 — > Xg-f) directly constrains \Cjj{mi,)\'^ + \Cj^{mb)\'^ at the leading order. Due to the strong enhancement factor 
mg/mi, associated with single S^^^^^^ insertion term in (^^^^^(to;,), 62^^^^^ (~ 10^^) are more severely constrained 
than S2:^^^'^\ However, if the left-right mixing of scalar bottom quark 3^;^ is large 0.5), 62^'"^^^ is constrained to 
be order of 10^^ since the double insertion term 82^^^^'^ 5^^^^^*^ is also enhanced by rrig/mi,. The branching ratio 
Bg n+fi- in MSSM is given as 

Br(i3,^^V-) - ^f^miTBji|A,pyT3I^[(l_4m2)|CQ,(m,)-C^^(m,)p-F 

|CQ,(TOf,) - C'Q^imi,) + 2m{Cio{mi,) - C[„{mb))\''] (37) 

where m = m^/ms,. In the moderate and large tan/3 case the term proportional to (Cm — C(q) in Eq. (|37|l can be 
neglected. The new CDF experimental upper bound of Br(i?s /x+/i^) is 1.5 x 10~^ |4flj at 90% confidence level. 
We have the constraint 

^ICqAmw) - C'Q^imwW + ICQ^mw) - C'Q^{mwW ^ 1-2 (38) 

Because the bound constrains \Cq^ — Cq,\ {i—l, 2), ^ we can have values of |CqJ and \Cq.\ larger than those in 

constrained MSSM (CMSSM) with universal boundary conditions at the high scale and scenarios of the extended 

I — I 

minimal flavor violation in MSSM 23] in which \Cq,\ is much smaller than |CqJ. Just like the constraint from 
Br(i3 -''^57), 62^'"^^^ is also constrained to be order of 10^^ by Br(i3s /i^/i^), if (Jgg^ is order of 0.5. At the same 
time we require that predicted Br{B — > Xgii^/j,^) falls within 1 a experimental bounds, which gives no new limits on 
parameters once the updated CDF bound of Br(i?.s /i+/i^) is imposed. It is shown recently that with the old CDF 
Br(_Bs — > ^'^ upper bound, 2.6 x 10^^ 42], the present experimental limit Rk [Rk = Br{B K fi^ fi^)/ Br{B — > 
Ke^e~)) < 1.2 puts constraints on Cq^ ^ which are similar to ones from Br(_Bs — > /I'^fi^) l32| and Higgs penguin 



Isi a 

is 0,1 



contributions (i.e., the terms relevant to Cq| ^) to Br(_B Xgfi^ ^, ) is order of 10% or less laiQ We obtained 
smaller contributions by calculations with the updated CDF bound. 

We also impose the current experimental lower bound AAfg > 14.4ps~^ ji^. The correlation between S^k and 
AAfs has been extensively discussed in the literature, in particular, in the fourth paper of ref.|23|. So in this paper 



Cq' ^ are the Wilson coefficients of the operators Qi'2 which are Higgs penguin induced in leptonic and semileptonic B decays and their 

definition can be found in Ref. |44| . By substituting the quark-Higgs vertex for the lepton-Higgs vertex, it is straightforward to obtain 
Wilson coefficients relevant to hadronic B decays. 



we just analyze the constraints on parameters from the lower bound. Because 62^^^^^ is constrained to be order of 
10"^ by Br{B — > Xgj), their contribution to AM^ is small. The dominant contribution to AMs comes from 62^'"^^^ 
insertion with both constructive and destructive effects compared with the SM contribution. Too large a destructive 
effect is ruled out, because SM prediction, AM^*^ — 17.3lto'7 13' '^^^Y slightly above the present experiment 

LL( RFC] 9 

lower bound. However 822, are constrained to be order of 10~ by the combined experimental measurement of 
Br(i? — > Xs^) and upper bound of Br(i?s — > fi'^n'). Their effects to AMg are limited. And we have checked, the 
effects are negligibly small with only one kind of chirality, LL or RR, while AMg can be enhanced to 25ps~-'^ with 
both kinds of chirality, LL and RR, however it is not strongly correlated with Si^Ks and Sr/'Ks provided that the 
other experimental constraints, in particular, those from Br(i3 —>■ X^"/) and upper bound of Br(_Bs — > have 
been imposed. 

As pointed out in Sec. II, due to the gluino-sbottom loop diagram contribution and the mixing of NHB induced 
operators with the chromomagnetic dipole operator, the Wilson coefficients Cgg can be large, which might lead to a too 



large Br(B Xsg). So we need to impose the constraint from experimental upper bound Br{B Xgg) < 9% 146 
A numerical analysis for Cg^^O has been performed in Ref.js^l- We carry out a similar analysis by setting both Csg 
and Cgg non-zero. 



B. Numerical results 

In the numerical calculations, we employ the latest Light-Cone Sum Rules results [sj] for the form factors of 
B K* , other parameters can be found in ref. [2^ . 

Before moving to numerical results, we discuss some unique features of i? ^ VV process. The contributions of 
non-primed operators to the helicity amplitude iJ+ are much smaller than those to while the contributions of 
primed operators to the helicity amplitude H- are much smaller than those to H^, because of the helicity flip of 
quarks and anti-quarks coming from non-primed or primed operators when they consist of a vector meson with some 
definite helicity. That is, in the transverse basis, Aq and are proportional to C — C", while Aj_ is proportional to 
C + C. Therefore, we have ~ \{C - C')/ {C + C')\. 

In numerical analysis we fix nig = niq = 500GeV, tan/5 = 10 and S^^^ = 0.4. We vary the NHB masses in the 
ranges of 91GeV < mh < 135GeV,91GeV < niH < 200GeV with ruh < mn and 200GeV < tua < 250GeV for the 
fixed mixing angle a = 0.6,7r/2 of the CP even NHBs and scan 623^ in the range < 0.06 for A=B and 0.01 

for ^ 7^ B (A = L, R). 

Numerical results for the correlation between longitudinal polarization /l and branching ratio Br(_B — > 4>K*) 
are shown in Fig.s QJEl where the correlation between oi B (/)K* and the indirect CP asymmetry Sb^4,k 
is also given. Fig. [H El El E are the results of insertions of (5^f ^, both 5^^^ and 5^^^ , (5^f ^, both 5%^ and 
62^^, respectively. In all four cases, /l can be dragged as low as 0.5, but the Br(i3 — > (j>K*) is smaller than the 
experimental measurement when ^ 0.5. On the other hand, there are some parameter regions with as low as 
0.5 and Sb^c/^k near 0.4, which is consistent with the present experimental measurements. In the case of new physics 
contributions from LR, RL insertions as shown in Fig. ^ El the only large effects come from the SUSY contributions 
of the chromo-magnetic dipole operator Qsg (and/or Q'^g) since the Wilson coefficient Cg™(TOf,) can be significantly 
larger than C|^(mf,). Because Qsg does not contribute to h = ±1 amplitude, only the longitudinal amplitude can be 
largely modified. The experimental measurement of ^ 0.5 requires that the magnitude of longitudinal amplitude 
in MSSM must be smaller than that in SM. Therefore, Br(_B 4>K*) in MSSM decreases, compared with SM, when 
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: 0,4 0.6 O.S 



FIG. 1: The correlations between /l and Br(B — > 4>K*), Sb—><pk with ^ll*^ insertion. The Br is in unit of 10 
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FIG. 2: The correlations between /l and Br(B ct>K*), Sb^-pk with (Jfl"^ and (Jfl*-^ insertions. 



/l ~ 0.5, as can be seen from Fig.^ [S] 

In the case of LL, RR insertions, the Wilson coefficient C'g^'(TOvi') could be largely modified and c['^{7) and c['^{7) 

could be large. Running from the mw scale, Cis^C'i^) can induce sizable C'sg (C'sg) '^fc scale, which effect we 

have discussed above. Running from a large electro- weak scale to mi,,Ci3{C'i^) can also induce large Ci3_i6 {Ci^^ig)- 

However, the updated CDF bound of Bs — > ^'''M^ imposed a stringent constraint on C'li which leads to 

-.(') 



that the Wilson coefficients of NHB induced operators are small and Cg^ can be largely modified only for small and 
moderate tan/3. The penguin insertions of operators QiI_iq as well as Qil 12 have been calculated and given in 
eqs.(13) and (14), but numerically Qi3_ig contributions to the magnitude of h — ±1 amplitude are small compared 
with the magnitude of /i = amplitude, due to small GF^{1) function in eq.(14). QiI^iq can contribute through 
tree-level to ai4„i6 in eq.(14) and Qil 12 also do. However, c[i are not large enough to enhance the transverse 
amplitudes sizably due to the constraint from Bs — » fi^ 11^ , as pointed above. So even though A±_ has the structure 
of C + C", which is different from the Ao.|| amplitudes, it is still impossible to fine-tune the magnitude of Aj_ to the 
level of \Aq\. Therefore, the overall contribution of LL, RR insertions are very similar to LR, RL insertions, as we see 
from Fig.sQ}^ 



The numerical results are obtained for rria 



-500 GeV. For smaller gluino and squark masses, the Wilson 



coefficient Cg^ becomes larger, which could have larger effect on the b to s transitions. However, the effect is indeed 
limited due to the constraint from B Xgg. For fixed rrig, the Wilson coefficient Cg^ is not sensitive to the variation 
of the mass of squark in the range about from 100 GeV to 1.5 TeV. Therefore, the numerical results are not sensitive 
to the squark mass and would have a sizable change when the gluino mass decreases. When the gluino and squark 
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FIG. 3: The correlations between /l and Br(i3 (pK*), Sb->4>k with 523*^ insertion. 



masses approach to infinity (indeed, the several TeV is big enough), SUSY effects drop, i.e., one reaches the decoupling 
limit. 

Before concluding, we will comment on two channels, B K*^ and B K*l~^l~, which share the same B K* 
form factors as i? ^ They have already been calculated within QCD factarization in ref. and new 

physics effects have been discussed in ref. [s^. For B — > K*j channel, the SM prediction is about 2 times larger 
than the experimental measurement. A way out to reduce the theoretical prediction of Br(i3 K*^) is to decrease 
the transverse form factors associated with B K* . Then the magnitude of transverse amplitude of B — > K*(l) will 
be decreased as well, and the polarization problem becomes even worse within the SM. We carry out an analysis of 
the correlations between Br(i3 K*j) and the polarization oi B ^ K*<j) within the new physics framework as we 
discussed above. We find that both Br(i? K*^) and can be accommodated within la limits only in the case 
of both LR and RL insertions as shown in Fig. However, in all the cases, the predicted Br(i3 K*(p) is still 
small when approaches 0.5. This situation can be relaxed to some extent in all the cases of insertions when we 
consider the B — > K* form factors ^|| and as defined in 0|, with 50% uncertainties. As an example, our results 
of the correlations between and Br{B — > K*(j)) are given in Fig. [Sja in the case of both LL and RR insertions. 
At the same time, Bi{B K*^) and Jl can be accommodated within Ict limits in all the cases of insertions. The 
situation oi B ^ K*l^l^ is more inconclusive due to the branch ing ratio measurement by BaBar and Belle with large 
uncertainties, and theoretically it has been discussed in ref. [siils^. 



V. CONCLUSIONS AND DISCUSSIONS 



In summary we have analyzed the B — > (pK* polarization puzzle in MSSM. The hadronic matrix elements of the 
new operators in MSSM for the decays have been calculated in the QCDF approach up to the order. Using the 
Wilson coefficients in ref. and hadronic matrix elements obtained, we have calculated the polarization fractions 
and branching ratios for the decays B — > (j)K* . It is shown that in the reasonable region of parameter space where 
the constraints from Bg — Bg mixing , B As7, B Xsg, B Xsji^ and Bg fJ-^l^^ are satisfied, the 
polarization fractions of the decays can agree with experimental data within la deviation. In particular, the puzzle 
for polarization in B ^ (pK* can be explained, while not in contradiction to the measurements of other two vector 
final states, in a large region of parameter space because we have included the contributions of the primed operators 
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FIG. 4: The correlations between /l and Br(B <t>K*), Sb->4,k with both 823^ and 523*^ insertions. 








FIG. 5: (a) The correlations between /l and Br(i3 K*^) with both 523'^ and 523*^ insertions, (b) the correlations between 
/l and Br(B K''4>) with both S^a^ and S^^^ insertions, while the B K* form factors and with 50% uncertainties. 



and new operators including the as corrections of hadronic matrix elements of them in MSSM. However the branching 
ratio is smaller than the measurements when the longitudinal fraction is near 0.5. We may not worry about it too 
much at present due to the large uncertainty in calculating hadronic matrix elements of operators. 

It is necessary to make a theoretical prediction in SM as precise as we can in order to give a firm ground for 
signalin g n ew physics. The twist-3 and weak annihilation contributions to i? ^ 0^^s in SM have been calculated 
in Ref. [2^ using the method in Ref. by which there is not any phenomenological parameter introduced. The 
numerical results show that the annihilation contributions to the decay rates are negligible, the twist-3 contributions 
are also very small, smaller than one percent. We expect that the conclusion would qualitatively remain for B — > (j)K* 
in MSSM, so that we have neglected the annihilation contributions in numerical calculations. 

In conclusion, we have shown that the recent experimental measurements on the polarization fractions in i? ^ (f'K* y 
which is difficult to be explained in SM, can be explained in MSSM if there are flavor non-diagonal squark mass 
matrix elements of second and third generations whose size satisfies all relevant constraints from known experiments 
{B Xsj, Bs — > /x^/i^, B — !■ XsfJ-^ fJ.^ , B — !■ Xgg, AMs, etc.). Therefore, if the present polarization puzzle persists in 
the future, it will be a signal for new physics beyond the SM and MSSM will be a possible candidate of new physics. 
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Appendix. Vertex and hard scattering contributions 

In the Appendix we give the explicit expressions of vertex corrections and hard scattering contributions at the as 
order for NHB induced operators which are not given in the content. 

The hard spectator contributions up to the leading twist in Eq. H14|l are given as following: 

^1344,15 = 



Tj± 4:Tr ^ JB J4, Jk* f df f du ^ f dv 

Nc -^T(l±75) Jo ^ Jo U ^ Jo V 

where is defined in ref. ^34^ and normalized as 







du (f>±{u) — 1. 

The vertex corrections up to the leading twist in Eq. (|14|l are as follows. 



''13,14,15 — ^ 



V^, = [121n h / du g{u) (f>±(u)] + the regularization scheme dependent constant, 

6 Jo ' 

o /I - 2a; , 

gix) — 3 mx — ITT 

\1 — X 

where we have used that (j)±{u,u) is symmetric with respect to u, u. Omitting regularization scheme dependent 
constants, we have 

= ^it = -12V^M- 



of a,. That is, -r^^ts = 0{al) 



We have verified that the dependance of = — |af4 + afs + \o.^q in eq.(ll) has been cancelled up to the order 

References 



[1] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 91, (2003) 171802 |hep-ex/0307026] . 
[2] Aubert et al. (BABAR Collaboration), Phys. Rev. Lett. 93 (2004) 231804 |hep-ex/04080T7| . 
[3] The Belle Collaboration, K. F. Chen et al, Phys. Rev. Lett. 94 (2005) 221804 |hep-ex/05030T3| . 
[4] B. Aubert et al. [BABAR Collaboration], hep-ex/0408093 

[5] K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 95 (2005) 141801 |hep-ex/0408T02l . 
[6] A. L. Kagan, Phys. Lett. B601 (2004) 151 .hep-ph/0405134. . 



14 



[7] M. Beneke et al., Phys. Rev. Lett. 83(1999) 1914; Nucl. Phys. B591(2000) 313. 

[8] H.-n. Li and H. Yu, Phys. Rev. Lett. 74 (1995) 4388; H.-n. Li and T. Yeh, Phys. Rev. D56 (1997) 1615; Y. Y. Keum, 
H.-n. Li and A. I. Sanda, Phys. Lett. B504 (2001) 6; Phys. Rev. D63 (2001) 054008. 

[9] H.-n. Li and S. Mishima, hep-ph/0411146 
[10] P. Colangelo, F.De Fazio, T.N. Pham, Phys. Lett. B 597, 291 (2004). 
[11] H.Y. Cheng, C.K. Chua, A. Soni, Phys. Rev. D71 (2005) 014030 |hep-ph/04093T7l ]. 
[12] H.-n. Li, |hep-ph/0411305| 
[13] A. L. kagan, |hepiph/0407076| 
[14] E. Alvarez et. al., |hep-ph/0410096| 
[15] A. .K. Giri and R. Mohanta, hep-ph/0 412107| 
[16] C. S. Kim, Y.-D. Yang, hep-ph/0412364" 
[17] H. -Y. Cheng, K. -C. Yang, Phys. Lett. B 511 (2001) 40. 
[18] P. K. Das, K.-C. Yang, Phys. Rev. D71 (2005) 094002 |hep^ph/0412313| . 
[19] W.-S. Hon, M. Nagashima, |hep-ph/0408007l 

[20] Y.-D. Yang, R.-M. Wang, G.-R. Lu, Phys. Rev. D72(2005) 015009 [hep^ph/041121l| . 
[21] C.-H. Chen, C.-Q. Geng, |hep-ph/0504145| 

[22] S. Baek, A. Datta, P. Hamel, O. F. Hernandez, D. London, Phys. Rev. D72 (2005) 094008 >hep-ph/0508149". 
[23] G.L. Kane, P. Ko, C. Kolda, J.-H. Park, H. Wang, L.-T.Wang, Phys. Rev. Lett. 90, 141803 (2003) ; Phys. Rev. D70 
(2004) 035015. 

[24] C.-S. Huang and Q.-S. Yan, Phys. Lett. B442 (1998) 209; C.-S. Huang, W. Liao and Q.-S. Yan, Phys. Rev. D59 (1999) 
011701; S.R. Choudhury and N. Gaur, Phys. Lett. B451 (1999) 86; K.S. Babu and C. Kolda, Phys. Rev. Lett. 84 
(2000) 228; C.-S. Huang et al., Phys. Rev. D63 (2001) 114021; ibid. 64 (2001) (2001) 059902(E); P.H. Chankowski, L. 
Slawianowska, P.R. D63 (2001) 054012; C. Bobeth et al., Phys. Rev. D64 (2001) 074014; hep-ph/0204225 G. Isidori, A. 
Retico, JHEP 11 (2001) 001; S. Baek, P. Ko and W. Y. Song, Phys. Rev. Lett. 89 (2002) 271801 ; JHEP 0303, 054 (2003). 
For a recent review, see, e.g., C. Kolda, *Tsukuba 2004, Supersymmetry and unification of fundamental interactions* 581, 
hep-ph/0409205 . 

[25] J.-F. Cheng, C.-S. Huang and X.-H. Wu, Phys. Lett. B 585 (2004) 287 hep-ph/0306086. . 
[26] J.-F. Cheng, C.-S. Huang and X.-H. Wu, Nucl. Phys. B701 (2004) 54. 
[27] C.-S. Huang and X.-H. Wu, Nucl. Phys. B657(2003) 304. 
[28] M. Beneke et al., Nucl. Phys. B6O6(2001) 245. 

[29] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68(1996) 1125 |arXiv:hep-ph/9512380| . 

[30] J.A. Bagger, K.T. Matchev and R.J. Zhang, Phys. Lett. B412(1997) 77; M. Ciuchini et al., Nucl. Phys. B523(1998) 501; 

C.-S. Huang and Q.-S. Yan, hep-ph/9906493 A.J. Buras, M. Misiak and J. Urban, Nucl.Phys. B586 (2000) 397. 
[31] F. Borzumati, C. Greub, T. Hurth and D. Wyler, Phys. Rev. D62(2000) 075005 . 
[32] G. HiUer, F. Kriiger, hep-ph/0310219 

[33] S. Wandzura and F. Wilczek, Phys. Lett. B 72 (1977) 195. 

[34] P. Ball, V.M. Broaun, Phys. Rev. D58, 094016 (1998); P. BaU, R. Zwicky, Phys. Rev. D71, 014029 (2005). 
[35] I. Dunietz, H.R. Quinn, A. Snyder, W. Toki, H.J. Lipkin, Phys. Rev. D43, 2193 (1991). 
[36] B.V. Geshkenbein and M.V. Terentev, Phys. Lett. 117(1982) 243 . 

[37] R. Kaiser and H. Leutwyler, Proceedings of Workshop on Nonperturbative Methods in Quantum Field Theory, Adelaide, 
1998 hep-ph/9806336 ; 

R. Kaiser and H. Leutwyler, Eur. Phys. J. C 17 (2000) 623 'hep-ph/0007101|. 
[38] H.-Y. Cheng, Y.-Y. Keum and K.-C. Yang, Phys.Rev.D65 (2002) 094023. 
[39] M. Ciuchini et al, J. High Energy Phys. 07(2001) 13. 
[40] The CDF Collaboration, |hep-ex/0508036l 



15 



[41] P. H. Chankowski and L. Slawianowska, Eur. Phys. J. C33 (2004) 123. 

[42] CDF Collaboration, F. Abe et al., Phys. Rev. D57 (1998) 3811. 

[43] A. Ali, P. Ball, L. T. Handoko, G. Hiller, Phys. Rev. D61 (2000) 074024. 

[44] Y.-B. Dai, C.-S. Huang and H.-W. Huang, Phys. Lett. B390(1997) 257. 

[45] A. Stocchi, Nucl. Phys. Proc. Suppl. 117(2003) 145 arXiv:hep-ph/0211245 . 

[46] A. Kagan, hep-ph/9806266 T.E. Coan et al. (CLEO Collaboration), Phys. Rev. Lett. 80 (1998) 1150. 

[47] J.-F. Cheng and C.-S. Huang, Phys. Lett. B554(2003) 155 . 

[48] R. Harnik, D. T. Larson, H. Murayama, A. Pierce, |aTXiv:hep-ph/0212180| 

[49] T. Feldmann, P. KroU and B. Stech, Phys. Rev. D 58(1998) 114006 'he p-ph/9802409] ; 

T. Feldmann, Int. J. Mod. Phys. A 15(2000) 159 "hep-ph/9907491 . 

[50] S. W. Bosch and G. Buchalla, Nucl. Phys. B 621 (2002) 459; A. Ah, A. Ya. Parkhomenko, Eur. Phys. J. C 23 (2002) 89. 

[51] M. Beneke, Th. Feldmann, D. Seidel, Nucl. Phys. B 612 (2001) 25; Eur. Phys. J. C 41 (2005) 173. 

[52] F. Kruger and J. Matias, Phys. Rev. D71 (2005) 094009. 



